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Abstract This paper deals with the analysis of sunspot number time series using 
the Hurst exponent. We use the rescaled range (R/S) analysis to estimate the 
Hurst exponent for 259-year and 11360-year sunspot data. The results show a 
varying degree of persistence over shorter and longer time scales corresponding to 
distinct values of the Hurst exponent. We explain the presence of these multiple 
Hurst exponents by their resemblance to the deterministic chaotic attractors 
having multiple centers of rotation. 
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1. Introduction 

Strong magnetic field present in the sun's outer regions is manifested by com- 
plex temporal dynamics, e.g., sunspots, solar wind velocity and solar flares. 
Magnetic activity manifests itself most clearly in sunspots. It has been found 
that chromosphcric flares show a very close statistical relationship with sunspots 
(Bray and Loughhead, 1979). Long term variations of solar activity may cause 



climatic changes on earth, whereas short term variations may be accompanied by 
fluctuations of certain meteorological parameters ( jWittmann, 1978] ) . Besides its 
11-year fundamental periodicity, solar activity as measured by relative sunspot 
number shows quasipcriodic variations with period ranging from 2 to 1100 years 
dMichelson, 1913||Kmuira, 1913[|Turner, 1913[|Kiral, 1961[|Zhukov and Muzalevskii, 1969| | 



Cole, 1973D . 

Alongside several periodicities, the solar activity also exhibits irregular fluc- 
tuations and these fluctuations were first assumed to be determined by the short 
term variation with a random distribution (Ruzmaikin, Feynman, and Robins on, 1994[ ).[ 



The rediscovery of the grand minima of solar activity (Eddy, 19761 led to a re- 
examination of the nature of the non-periodic part of the variations of the sun's 
activity. Solar activity in the frequency range from 100 to 3000 years includes an 
important continuum component in addition the well-known periodic variations 



(Ruzmaikin, Feynman, and Robinson, 1994) 
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The Hurst exponent is a parameter that quantifies the persistent or anti- 
persistent (past trends tend to reverse in future) behavior of a time series. It 
determines whether the given time series is completely random or has some long- 
term memory. Ruzmaikin et al. (1994) examined whether or not the nonperiodic 
variations in solar activity are caused by a white-noise, random process. They 
evaluated the Hurst exponent for a time series of 14 C data from 6000 BC to 1950 
AD. They find a Hurst exponent of ps 0.8 indicating a high degree of persistence 
in the variations of solar activity. Kilcik et al. (2009) used the monthly ISSN 
(international sunspot numbers) data for the last 3000 years to evaluate the 
Hurst exponent with a view to predict the sunspot activity for solar cycle 24. 
Xapsos et al. (2009) used the reconstructed sunspot numbers for the past 11 360 
years by Solanki et al. (2004) to find the Hurst exponent of ~ 0.8 and also showed 
the evidence of 6000-year periodicity in the reconstructed sunspot numbers. In 
all the above studies involving the Hurst analysis to understand the persistent 
behavior of the sunspot data, a single Hurst exponent was estimated although 
there arc many scaling regimes which give different Hurst exponents. In this 
paper, we use the Hurst analysis on 259-year and 11 360-year data sets and find 
multiple Hurst exponents in each time series. We explain the presence of multiple 
Hurst exponents for a single time series using systems from deterministic chaotic 
dynamics with multiple centers of rotation. Our results conclude that estimating 
a single Hurst exponent from the data where different linear scaling regimes exist 
may be improper. 

In Section [2] we review the R/S method to calculate the Hurst exponent of a 
time series. The results for sunspot data are discussed in Section [3] Results for 
different chaotic models having one or more centers of rotation in phase space 
are given in Section 01 The conclusions are given in Section [5l 

2. Hurst Analysis and R/S Measure 

The R/S method to find the Hurst exponent was proposed by Mandelbrot and 
Wallis (|1969b ) which can be summarized as follows: 

Let Xi, i = 1,2, ...,7V, be an observed time series whose Hurst exponent is 
to be computed. Let us now choose a parameter w (temporal window such that 
Wt < w < N where wt is the Theiler window ( |Theiler, 1986[ )) and consider the 
subsets of the data Xi, i = to,t + l,t + 2, ...,to + w — 1, where 1 < t < 
N-w + 1. 

We then denote the average of these subsets as 



Let S(to,w) be the standard deviation of Xi, during the window w i.e., 




t a +w~ 1 



i=t a 




l 

2 
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Next, new variables j/j, i = 1, 2, w, and range R are defined as 



to+i-l 

y i (t Q ,w) = ^2 [xk - x(t ,w)], 

k—to 



R(t ,w)= max j/i(t ,w)- min yi(t ,w), 



which allows one to define Reseated range measure R/S as 

R(t ,w) 



(R/S)(t ,w) 



S(t ,w) 



Taking To = 1,2, N — + and computing (R/S)(to,w) for time lag w 
the rescaled range for the time lag w is finally written as the average of those 
values 

R/S= M 1 — 52(R/S)(to,w). (1) 
TV — w + 1 

to 

It has been observed that the rescaled range (R/S) over a time window of 
width w varies as a power law: 

(R/S) w = kw H , (2) 

where k is a constant and H is the Hurst exponent. To estimate the value of 
the Hurst exponent, R/S is plotted against w on log-log axes. The slope of the 
linear regression gives the value of the Hurst exponent. If the time series is purely 
random then the Hurst exponent (H) comes out to be 0.5. If H > 0.5, the time 
series covers more 'distance' than a random walk, and is a case of persistent 
motion, while if H < 0.5, the time series covers less 'distance' than a random 
walk it shows the anti-persistent behavior in time series. However, for periodic 
motions, the Hurst exponent is 1. 



3. Sunspot Data and R/S Analysis 

Sunspot number (SSN) is continuously changing in time and constitutes a time 
series. Figure HJa) shows the monthly-averaged sunspot numbers from Sunspot 
Index Data Center SIDC (Ihttp: / / www.sidc.be/sunspot-datal) from 1749 up to 
the present. The power spectrum of this data set is shown in Figure^ a). We 
further estimate the Hurst exponent for this time series using the R/S method. 
Figure G^a) shows the log-log plot for the 259-year data showing the presence 
of multiple Hurst exponents for this time series. The first bending is due to the 
prominent 11-year cycle and the estimated Hurst exponent 1.13 in the linear 
regime before this bending. The next prominent bending is at roughly 90 years 
showing the Wolf-Glcissberg cycle. The estimated Hurst exponent in the second 
linear regime between 11 years and 90 years is 0.83. 
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Figure 1. a) Monthly sunspot record from January 1749 to December 2008. b) Reconstructed 
sunspot record from Solanki et al. (2004). The time period covers 11 360 years with observations 
in 10-year increments. 
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Figure 2. Power spectra for a) monthly sunspot record from January 1749 to December 
2008. b) Reconstructed sunspot record from Solanki et al. (2004). The time period covers 
11 360 years with observations in 10-year increments. 



We further analyze the reconstructed sunspot time series of 11360 years 
(Solanki et al., 2004) shown in Figure [ljb). Figure EJb) shows the power spec- 
trum of this time series indicating the presence of two long-term periodicities at 
nearly 2300 and 6000 years. The R/S analysis also shows two bondings around 
these periods (Figure [3jt>) ) . Again, for this data set we get two linear regimes, 
the Hurst exponent of the first part being 0.77 and the second part being 1.26. 

If we were to calculate a single Hurst exponent for these two time series we will 
get H w 0.9 for 259-year and H « 0.8 for 11 360-year data (Figure^a) and (b) 
respectively). These values of the Hurst exponent agree with the previously esti- 
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Figure 3. log(R/S) vs. log(ui) for sunspot number a) from January 1749 to December 2008 
b) reconstructed sunspot record of 11 360 years. The Hurst exponents with corresponding 
standard deviations are a) 1.13±0.003; 0.83±0.003 b) 0.77±0.001; 1.26±0.004. 



mated values of Hin the literature (Ruzmaikin, Feynman, and Robinson, 1994; 



Gil-Alana , 2009 Xapsos and Burke, 2009). However, as we have demonstrated, 



the log-log plots of sunspot data show two remarkably distinct scaling regimes 
and hence estimating only one Hurst exponent may be improper. In the next 
section we examine in detail the reasons for these different scaling regimes by 
using examples of a variety of time series of standard chaotic systems. 



4. Chaotic Models and R/S Analysis 

In order to understand the behavior of different dynamical systems using the 
Hurst exponent, we take two extreme examples of time series generated from 
random and periodic motion respectively. First we take a random time series 
with 100 000 data points distributed uniformly in the range [0, 1]. Using the R/S 
method we compute the Hurst exponent as discussed in Section 2 and find the 
Hurst exponent H = 0.5 as given by the slope of the line in Figure H^a), which 



is expected for a purely random time series ( |Mandelbrot and Wallis, 1969a[ ). 

We further consider the other extreme case of a purely periodic signal of 
period one generated from sin(x), with x in the range [0,1000] with step size 
0.01. Figure [4jb) shows the plot of \og(R/S) vs. log(w) and the slope of the 
linear regime in this case is 1. One prominent difference between these two 
extreme cases of a purely random signal and a periodic signal (Figure H£a) an d 
(b) respectively), is that for a purely random signal, the plot of log(R/S) vs. 
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log(w) has a constant slope, while for a periodic signal it gets saturated and 
starts oscillating after a certain value of w (marked by an arrow). The value 
where the bending begins corresponds to the period of oscillation T = 2-7r/0.01, 
and has been verified by the power spectrum. The bending, therefore, gives us 
information about the frequency of the given signal. 



4.1. Chaotic Time Series 



In the past few decades, chaotic motion, observed in deterministic systems, has 
received great attention due to its presence in systems from physical, chemical, 
biological, ecological, physiological to social sciences. These chaotic motions are 
temporally aperiodic and are strange because they have a fractal geometry. Since 
chaotic motion is neither periodic nor random we can expect its behavior to be in- 
between two extreme cases of randomness and periodicity. Therefore, we can also 
expect the Hurst exponent to be different from 0.5 and 1. An example of a system 
exhibiting chaotic motion is the celebrated Rossler attractor ( jRossler, 1979[ ), 

x = -y- z, 
y = x + ay, 

z = b + z(x — c), (3) 

where a, b and c are control parameters. A chaotic trajectory (with a = 0.1, 
b = 0.1 and c = 18.0) in the x — y plane is shown in Figure (5ja) which rotates 
around an unstable period-one fixed point ((c — d)/2, (— c + d)/2a, (c — d)/2a) 
where d = Vc 2 — 4afe . The plot of log(R/S) vs. log(u>) (Figure E^c)) shows 
that there is linear regime (before marked arrow) after which it gets saturated. 
The slope of linear region gives H « 1 . The first bending of the above-mentioned 
curve gives the frequency which matches with the frequency obtained from power 
spectrum or peak to peak analysis of the amplitude. This system is dissipative 
and has to be bounded in the sub-phase space flRossler, 1979 ) and the bending 



here corresponds to the folding nature of the dynamics. 
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In order to see if the behavior is replicated in another similar system, we 
consider the Chua oscillator (Chua et at, 1993), 

x = ci(y- x-p(x)), 
i) = c 2 (x -y + z), 

z = -C3V, (4) 

where p(x) is defined as, 

p{x) = m\X + ((m — mi)(| x + 1 | — | x — 1 |)) /2. 

We take c\ = 15.6, C2 = l,mo = —8/7, mi = 5/7 and first select the parameter 
C3 = 33 such that the motion will be a single-scroll type as shown in Figure [5]^b) 
(there is a symmetric attractor also, depending upon the initial conditions). The 
log(R/S) vs. log(iu) plot for variable x is shown in Figure [5]Jd) and the Hurst 
exponent is found to be H as 1. This also shows saturation at the average time 
period (shown by arrow). These two examples of chaotic dynamics, Rossler and 
Chua systems, clearly demonstrate that whenever there is a single center of 
rotation the log(i£/£) vs. log(w) plot shows a linear scale only up to the average 
period of the attractor after which it saturates. 

In nature one may come across many dynamical systems which are chaotic 
and their trajectories rotate around more than one center of rotation. In order to 
see the behavior of log(R/S) vs. log(w) plot, we first consider the Lorenz system 
dLorenz, 1963] ), 

x = a(y-x), 
V = x(p - z)-y, 

i = xy — j3z. (5) 

Shown in Figure [6ja) is the chaotic trajectory in the x — y plane at parameter 
values a = 16, p = 50 and f3 = 4. This clearly shows that the trajectory rotates 
for some time around the fixed points: 

(V/3(P"1). vW- 1 ), (P - !)) and (-y/P(p-l), -V/3 (/>-!)* (P - !))> 

while some time it rotates around all of these fixed points (including (0,0,0)). 
The \og(R/S) vs. log (it;) plot in Figure [B](c) shows that there arc two regimes 
of linear scaling. The slope of the first part of the linear region gives H — 0.93 
and the second part gives H = 0.64. The first linear regime corresponds to the 
trajectory rotating about individual fixed points while the second is for all three 
fixed points. This indicates that the dynamics around individual fixed points is 
different from that of the combined one. Therefore taking the slope of individual 
regimes can give more details of the intrinsic dynamics. This method also allows 
us to estimate intrinsic frequencies of the system. 

In order to see this behavior in another system where chaotic motion contains 
many unstable fixed points around which the trajectory revolves, we consider 
the Chua circuit represented by Equation Figure E{b) shows the trajectory 
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Figure 5. Plot of trajectories in the x — y plane for chaotic motion of (a) Rossler, b) Chua 
oscillators having single center of rotation. Corresponding plots of log(R/S) vs. log(ui) are 
given in panels c) and d) respectively. 



having a double-scroll chaotic motion (Chua system at C3 = 28). Similar to the 
Lorenz system, this system confirms the existence of two different regimes of 
linear scaling having the Hurst exponents, H = 1 and 1.1. These two examples 
confirm that for multiple centers of rotation, there are many regimes of linear 
scaling for which distinct Hurst exponents can be estimated. 



5. Conclusions 

In this paper, we use the Hurst analysis on 259-year and 11 360 -year data sets 
and find multiple Hurst exponents in each time series. We explain the presence of 
multiple Hurst exponents in a single time series using systems from deterministic 
chaotic dynamics with a single center of rotation ( Rossler and single-scroll Chua 
oscillators) as well as multiple centers of rotation (Lorenz and double-scroll Chua 
oscillators). We have shown that in the sunspot data, two distinct linear scaling 
regimes exist for which two distinct Hurst exponents could be estimated implying 
a variety of persistent behavior. The results show very high degree of persistent 
behavior till 11 years (H w 1.1), a slightly less persistent behavior till 100 years 
{H w 0.8), and comparatively less persistent behavior till 2300 years (H ss 0.77). 
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Figure 6. Plot of trajectories in the x — y plane for chaotic motion of a) Lorenz and b) Chua 
oscillators having multiple centers of rotation. Corresponding plots of \og(R/ S) vs. log(to) are 
given in panels c) and d) respectively. 
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